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Abstract 

We present a calculation of the D* + , D + and D° meson single inclusive produc- 
tion cross section for the Tevatron Run II. We use the FONLL approach in pertur- 
bative QCD, which, besides including the known next-to-leading order corrections, 
also provides for the resummation at the next-to-leading logarithmic level of terms 
enhanced at large p T by powers of log(p T /m), where m is the charm mass and p T 
is its transverse momentum. Non-perturbative effects in charm hadronization are 
extracted, in moment space, from recent ALEPH data for D* fragmentation in 
e + e~ collisions. 



1 Introduction 



Heavy flavour production at hadron colliders has been mostly studied for bottomed 
hadrons so far, due to the possibility of using muons or J/ty to tag the B hadrons. At 
CDF, thanks to the development of the silicon vertex trigger, charm production studies 
have become possible. This opens up the interesting possibilities of comparing QCD 
predictions for charm production with experimental results, much in the same way as 
one does for bottomed hadrons. One is then especially interested in seeing whether the 
predicted charm cross section bears the same relationship with respect to data as the 
bottom cross section. 

In the present work, we provide cross sections for charm production for the Tevatron 
Run II. The calculation of the cross section is performed according to the FONLL 
prescription |L, and the effect of a non-perturbative fragmentation function is included 
in full analogy with the work of ref. [2|, where a prediction for B production at the 
Tevatron Run I is given. We will thus not review the FONLL framework, and the 
moment space method for the extraction of the non-perturbative fragmentation function, 
rather referring the reader to refs. PJ E] 1 The only novel aspect of the present work 
that needs to be discussed is the phenomenological extraction of the non-perturbative 
.D-meson fragmentation function, which will serve as input to our calculation of the 
D spectrum. Good data on D* fragmentation have been published by the ALEPH 
collaboration jH] , together with a study yielding the relative fraction of directly produced 
D and D* mesons. In section we will first describe how we have constructed the D° 
and D + fragmentation functions in terms of the D* one and of the relative fractions 
of directly produced vector and pseudoscalar states. In section |3] we discuss how the 
D* non-perturbative fragmentation function itself was extracted from the ALEPH data. 
In section 0] we present our phenomenological predictions, and in section |3] we give our 
conclusions. 

2 Construction of the D° and D + non-perturbative 
fragmentation functions 

Non-perturbative information concerning the hadronization of the charm quarks into 
the observed D mesons must be extracted from experimental data. High quality data at 
present only exist for fragmentation into the vector state D*, as measured by the ALEPH 
collaboration [3| at LEP. The extraction of this non-perturbative fragmentation function 
(FF from now on) will be described in the next section. In the present one we show how 
this FF can be suitably modified in order to describe the production of the pseudoscalar 
states D° and D + too. 

The D° FF has been postulated to be a linear combination of £) 's coming from the 

Suffice here to say that the perturbative FONLL calculation includes the full next-to- leading order 
calculation for heavy quark production (and therefore also the power suppressed m/pT terms, all 
important in the pr ~ rn region), plus the resummation to next-to-leading logarithmic accuracy of the 
log^y/m) terms, which are large where px S> tn. The moment space method, on the other hand, takes 
care of properly extracting from experimental data, and importing in the process at hand, the relevant 
non-perturbative information. 
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decays of the D* 0, s and _D* +, s, and D 0, s not coming from D* decays, that we call here 
"primary" D 0, s : 

F(c -> D°) = F p (c -> Z?°) + F(c -> £>*+) <g> F(L>* + -> L>°) 
+ F(c^ D*°)®F(D* + ^ D°) 

where F p (c — > D°) is the fragmentation function for primary D° production, F(c — > D*°) 
and F(c — > are the fragmentation functions for D* states production, F(D*° — ► 

D°) and F(L>* + -> D°) describe the decay of a L>* into a D°. 

We assume that the production of primary D and D* is flavour independent, and the 
only flavour dependencies arise in the decay of the D*. In other words, one assumes as 
usual that isospin is a good symmetry, except that, due to the known fact that the 
decay mode is barely accessible to the D*, small isospin violations are amplified in the 
decay. Under this assumption, we can extract from the ALEPH paper |Hj the primary 
D° fraction and the D* + and D*° fraction 2 : 

BR(c^ D° p ) = 0.168, BR(c -> D*) = 0.233 . (1) 

Observe that these values conflict with the naive assumption of spin counting, that 
would predict a ratio of the D to D* states equal to the ratio of the number of spin 
states, 1/3. 

We now need a model for the D* FF. For sake of definiteness, we have taken the 
analytic forms of ref. 4J for the fragmentation into pseudoscalar (P) and vector (V) 
states. Such FF's will be denoted by D^ FY {z;r) and D^ FY {z;r) respectively, and 
understood to be normalized to one. A single non-perturbative parameter r describes 
both pseudoscalar and vector production in this formulation, and we can thus use the 
parameter fitted to D* production as input in the pseudoscalar (D° and D + ) fragmen- 
tation functions. We will comment in the following about the influence of this choice on 
our result. 

The D* decay into a D° could be calculated exactly using the known decay chains. 
However, since the decay D* — > ttD° is dominant and the decay products are nearly at 
rest in the decay frame, the effect of this decay is simply modeled as a rescaling of the 
longitudinal momentum by a factor m^/mu* in the FF for D* production: 

Dbcfy( z '-> r ) 



O(x) being the Heaviside step function. We can therefore write 

D^ D \z-r) = BR(c^D p )D^ FY (z- ir ) + 

+ BR(c^ D*)(BR(D*° -> D°) + BR(ZT+ -> D°))D^ FY (z;r) (3) 

2 Here and in the following, the branching fractions are presented and used without an associated 
experimental uncertainty. While such errors are in principle available from the experimental papers, 
we have refrained from quoting them as their contribution to the total uncertainty will turn out to be 
negligible. 
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where |S] 

BR(D*° -> D°) = 1, BR(ZT+ -> D°) = 0.677 . (4) 
Putting in numbers, eq. © reads 

D^ D \z-r) = VAQ%D^ FY (z-r) + Q.WD { ^ FY (z-r) (5) 

showing that about 2/3 of the fragmentation function depends in a precise way on the 
D* fragmentation. 

The FF for D + production can be constructed along similar lines. Since D*° cannot 
decay to D + we only have 

F(c D+) = F p (c -> D + ) + F(c -> £>*+) ® F(£>* + -> £)+) (6) 

and therefore 

D c -* D+ (z;r) = BR(c^D;)D^ FY (z;r) + 

+ BR(c^ £)*)BR(D*+ -> D + )D^ FY (z;r) (7) 

with 3 0E] 

BR(c -> D+) = 0.162, BR(D* + -> D + ) = 0.323 . (8) 

Numerically 

D c ^ D+ {z; r) = 0.162D { ^ FY (z; r) + 0.07153Dg8 fY («; r) . (9) 

For the sake of completeness, we now write down the FF for production of D* states: 

D^ D \z-r) = BR(c - D*)D^ FY (z;r) = 0.233D^ FY (z;r) . (10) 

The non-perturbative parameter r, which in the model of ref. j3] can be interpreted as 
the ratio of the constituent mass of the light quark to the mass of the meson, will be 
more precisely determined in the next section by a comparison with the ALEPH D* 
data. 

It is worth noting now that, in the present work, the model of ref. [1] is needed in 
order to compensate for the lack of an explicit measurement of the direct component of 
the fragmentation function. In the case of the D° production the influence of the 
model is marginal, since most D° are produced through D* decays. In the case of D + 
production we have the opposite situation. Therefore, our D + spectrum will strongly 
reflect this model dependence. New data on D + fragmentation may eventually help in 
reducing this uncertainty in the future. 

3 Extraction of the D* non-perturbative fragmenta- 
tion function 

The non-perturbative fragmentation function has been extracted from D* fragmentation 
data from the ALEPH collaboration We assume the measured energy fraction (x F ) 

3 Note that BR(c -> D+) = 0.162 is actually to be considered equal to BR(c -> D°) = 0.168 within 
the experimental uncertainties and consistently with our flavour-independence hypothesis. 
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Figure 1: a) Left plot: px spectrum of the charm cross section in pp collisions at y/s = 1960 
GeV, compared to A/pj, behaviours, b) Right plot: Mellin moments of D* fragmentation 
distribution as calculated from ALEPH data [3], compared to pure perturbative QCD and to 
pQCD convoluted with the non-perturbative form D^q FY N (r = 0.1). 



distribution (normalized to one) to be described by a convolution (a product, in Mellin 
moments space) of a perturbative contribution (which describes the production of a 
c quark in the hard interaction) and the BCFY non-perturbative form describing its 
hadronization into the D*\ 

a N (c - D*)(r) = Df CD D^ FYN (r) . (11) 

The non-perturbative parameter r is then determined by comparing <jjv to the measured 
( x e~ 1 ) values. 

According to the discussion published in [2] about B meson production at CDF, 
we have taken care that the fitted non-perturbative fragmentation function properly 
describes the Mellin moments of the D* spectrum in e + e~ experimental data around 
N = 4,5 rather than the whole xe shape, as this information is the most relevant 
one when calculating the hadronic production. Indeed, as shown in figure G^a), m the 
region of interest between — 5 and pr = 20 GeV, the charm transverse momentum 
distribution predicted by the FONLL calculation falls off like A/plp, with n between 4 
and 5. 

Figure db) shows a comparison of the Mellin moments calculated from the D* 
fragmentation data of the ALEPH Collaboration |3] with the result of a convolution 
of a perturbative and a non-perturbative distribution. The perturbative component 
jjpQ cd resum s to all orders and to next-to-leading logarithmic [H] accuracy the large 
q" \og k (Q 2 /m 2 ) (k < n) terms (Q being the e + e~ centre-of-mass energy, and m the 
heavy quark mass) which appear in the fixed order calculation. This resummation is 
analogous to the one performed in the FONLL result, where a" \og k (p^/m 2 ) terms are 
resummed. For the non-perturbative distribution, the functional form for vector states 
given by Braaten et al. [Ij has been chosen. Though no rigorous fit 4 of its free parameter 

4 The choice of not performing a rigorous fit of the moments is related to the fact that their errors have 
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Figure 2: QCD predictions for the D° differential cross section at the Tevatron Run II. 
CDF preliminary data are also shown for comparison. 



r has been performed, one can see from figure ^b) that the choice r = 0.1 gives a good 
description of the relevant moments around N = 4,5. The failure of the curve to describe 
higher- N moments is of course related to unaccounted-for sub-leading logarithms related 
to soft-gluon emission becoming more important. Resummation of these terms would 
improve the picture but, though available 0|S], it has not been included since it would 
not affect this phenomenological application 5 . 

4 Phenomenological results 

We show in figs. I2I3I4I our predictions for the da/dp T (\y\ < 1) cross sections for D°, 
D* + and D + production at the Tevatron Run II, obtained with the FONLL perturbative 

been only naively estimated from the errors of each xe bin, a full correlation matrix being unavailable. 
Lacking a rigorous calculation for the errors on the moments and the correlations between the various 
moments, a fit would yield results of little significance as far as the errors on the fitted parameters are 
concerned. 

5 By tweaking the factorization scales which appear in the resummed perturbative calculation one 
could improve the behaviour of the curve in the N > 6 region even without resumming the soft- 
gluon logarithms. This implies a different value for the non-perturbative parameter r when fitting the 
experimental moments. We have checked that the differences compensate, and the predictions for D 
meson production in pp collisions remain identical within a few percent. 
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Figure 3: Same as fig. El but for D* + production. 



result combined with the non-perturbative fragmentation functions given in eqs. (fiTty .flHj) 
and © respectively. The perturbative calculation is consistent with the one employed in 
calculating e + e~ fragmentation when fitting the D* non-perturbative FF, and it is also 
shown for comparison. The three curves in each plot are the lower, central and upper 
values of an uncertainty band produced by varying independently the factorization and 
renormalization scales in the transverse mass region m T /2 and 2m T . The charm mass 
has been fixed here to m = 1.5 GeV, and the CTEQ6M parton distribution functions 
set, corresponding to A^ 5 ' = 0.226 GeV, has been used. We have verified that accounting 
also for the effects of varying the QCD scale in the range 0.192 < A^ 5 ^ < 0.241 GeV and 
the charm mass in the range 1.2 < m < 1.8 GeV does not produce a larger band than 
the one given by factorization and renormalization scale variations 6 . 

The preliminary data measured by the CDF Collaboration |J| are also presented 
in the plots, superimposed to the theoretical curves. They can be seen to be in fair 
agreement with the predictions. The slopes in particular are very well described, and 
they differ from that of the QCD prediction for charm quarks, pointing to a correct 

6 An explanation of the limited sensitivity to the value of m is that the effects of mass variations are 
largely compensated by the need to also fit a different non-perturbative parameter r from the e + e~ 
ALEPH D* data, as m also enters the evolution of the perturbative calculation. It should also be 
stressed that the variation of either the mass or the QCD scale is not computed fully consistently, 
because we are not in a position to determine the effect of these variations on the parton densities. We 
therefore content ourselves with computing them at fixed parton densities. 
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Figure 4: Same as fig. El but for D + production. 



determination of the softening effects related to the non-perturbative fragmentation 
functions. Explicit numerical values for the theoretical predictions are given in Table EJ 
and the Data/Theory ratios are shown in fig. 0] together with the theoretical uncertainty 
bands, here defined as upper/lower theoretical prediction divided by the central value. 

The bands shown in the plots of figs. 121-3141 do not account for all the possible un- 
certainties. Besides the effects of renormalization/factorization scales, mass and QCD 
scale variations, previously discussed, one could also consider changing the Parton Dis- 
tribution Function set and/or modifying the model describing the hadronization of the 
charm quark into the mesons. At the present level of accuracy, however, we do not 
expect these issues to change significantly the picture presented in the figures. 

5 Conclusions 

We have presented predictions for D* + , D + and D° transverse momentum distributions 
as measured by the CDF Collaboration at the Fermilab Tevatron Run II in pp collisions 
at a centre-of-mass energy of 1960 GeV. 

The calculation consists of next-to-leading log resummed QCD (matched to a fixed 
next-to-leading order result) combined with non-perturbative information describing the 
charm hadronization into the mesons, and determined from moment-space analysis of 
D* + data measured by the ALEPH Collaboration in e + e~ collisions. 
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Vt (GeV) 


da/dp T (\y\ < 1) (nb/GeV) 




D° 


D*+ 




5.75 


.3851;^ xlO 4 


-i92±:5S xio 4 


.1701;^ xlO 4 




7837±220±884 


_ 


_ 


6.5 


.232^ xlO 4 


.iiTt:ffiS xio 4 


.1031;^ xlO 4 




4056±93±441 


2421±108±424 


1961±69±332 


7.5 


.124t;££ xlO 4 


•634i;^5 xlO 3 


.554™ xlO 3 




2052±58±227 


1147±48±145 


986±28±156 


9 


•539l; 24 | xlO 3 


.279i;^ xlO 3 


.2421;^ xlO 3 




890±25±107 


427±16±54 


375±9±62 


11 


•203±;^ xlO 3 


.107™ xlO 3 


.QlTt^S xlO 2 




327±15±41 


148±8±18 


136±4±24 


16 


•2921;^ xlO 2 


•157^035 x 10 2 


.133±:$g xlO 2 




39.9±2.3±5.3 


23.8±1.3±3.2 


19.0±0.6±3.2 



Table 1: Theoretical predictions, with uncertainties, for the differential cross sections 
of D°, D* + and -D + meson production in pp collisions at the Tevatron Run II. The 
preliminary experimental measurements from the CDF Collaboration J[| are also shown 
in the second lines at each p T value. The statistical and systematic errors have been 
added in quadrature in the figures. 



The overall agreement is good, with the data lying on the upper limit of the theoret- 
ical uncertainty band. This behaviour is similar to what observed in B mesons E] 
and 6-jets [TTJ [T2] production at the Tevatron. This consistency is encouraging, and 
many improvements, like tuning of scales and/or quark mass, adding higher order per- 
turbative contributions, or a more accurate determination of non-perturbative effects, 
could easily conspire to bring the theoretical prediction in even better agreement with 
the experimental data. It has been speculated that new physics mechanisms may be 
needed in order to explain the measured excess in the b cross section [T^j. The pres- 
ence of 'new physics' effects cannot certainly be ruled out at this stage. However, the 
fact that both charm and bottom production data bear similar relationship to the the- 
oretical prediction gives an indication that uncalculated QCD effects are a more likely 
explanation of the discrepancy. 
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